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Abstract. 

Interferometry with ultracold atoms promises the possibility of ultraprecise and 
ultrasensitive measurements in many elds of physics, and is the basis of our most 
precise atomic clocks. Key to a high sensitivity is the possibility to achieve long 
measurement times and precise readout. Ultra cold atoms can be precisely manipulated 
at the quantum level, held for very long times in traps, and would therefore be an 
ideal setting for interferometry. In this paper we discuss how the non-linearities from 
atom-atom interactions on one hand allow to efficiently produce squeezed states for 
enhanced readout, but on the other hand result in phase diffusion which limits the 
phase accumulation time. We find that low dimensional geometries are favorable, 
with two-dimensional (2D) settings giving the smallest contribution of phase diffusion 
caused by atom- atom interactions. Even for time sequences generated by optimal 
control the achievable minimal detectable interaction energy AE^^^ is on the order 
of 10~^/i, where /i is the chemical potential of the BEG in the trap. From there we 
have to conclude that for more precise measurements with atom interferometers more 
sophisticated strategies, or turning off the interaction induced dephasing during the 
phase accumulation stage, will be necessary. 
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1. Introduction 

Interferometry is the method of choice in achieving the most precise measurements, or 
when trying to detect the most feeble effects or signals. Interferometry with matter 
waves [Ij became a very versatile tool with many applications ranging from precision 
experiments to fundamental studies. 

In an interferometer, an incoming ^beam^ (light or matter ensemble) is split into 
two parts (pathways), which can be separated in either internal state space or real 
space. The splitting process prepares the two paths with a well-defined relative phase 
S9{t = 0) = 9i{t = 0) — 92{t = 0). After the splitting they evolve separately, and 
can accumulate different phases Oi{t) and 02{t) due to the different physical settings 
they evolve in. Finally, in the recombination process after time T the relative phase 
S9{T) = 9i{T) — 92{T) accumulated in the two paths can be read out. 

The sensitivity of an interferometer measurement depends now on two distinct 
points: How good can the phase difference 59 be measured, and how long can 
one accumulate a phase difference in the split paths. For perfect read-out contrast 
and standard (binomial) splitting and recombination procedures, the uncertainty in 
determining /\9 is given by the standard quantum limit = l/^/N^ where is 
the number of registered counts (e.g. atom detections). The second point concerns 
the question of how long the beams can be kept in the 'interaction region' of the 
distinguishable interferometer arms, and for how long the paths stay coherent. Ultra 
cold (degenerate) atoms can be held and manipulated in well controlled traps and 
guides, and therefore promise ultimate precision and sensitivity for interferometry. 
Both dipole traps [2j and atom chips (31 IH El E] have been used to analyze different 
interferometer geometries [Zl El El HD] and employed for experimental demonstration of 
splitting dH [121 [El [H] and interference [l5l[l6l[I7l[IHl[l9l[20l[ai[2a with 
trapped or guided ultra cold atoms. 

The power of interferometry lies in the precision and robustness of the phase 
evolution, which provides the measurement stick. This robustness of the phase evolution 
is based on the linearity of time propagation in the different paths, which is the case 
for most light interferometers, atom interferometers [Ij with weak and dilute beams, 
or neutron interferometers [25]. Measurements loose precision when this robustness of 
the phase evolution cannot be guaranteed. This is the case if the phase evolution in 
the paths depends on the intensity (density), that is, when the time evolution becomes 
non-linear. Atom optics is fundamentally non-linear, the non-linearity being created 
by the interaction between atoms. For ultra cold (degenerate) trapped Bose gases 
this mean field energy associated with the atom-atom interaction can even dominate 
the time evolution. Consequently, in many interferometer experiments with trapped 
atoms the atom- atom interaction creates a non-linearity in the time propagation, 
and the accumulated phase depends on the local atomic densities. Thus, number 
fluctuations induced by the splitting cause phase diffusion [26l E7j, which currently 
limits the coherence and sensitivity of interferometers with trapped atoms much more 



CONTENTS 



4 



then decoherence coming from other sources hke the surface [28l [29l l30] . 

In the present manuscript we wiU discuss the physics that leads to degradation of 
performance of an atom interferometer with trapped atoms, and how one can counteract 
it by using optimal input states. We first discuss the performance of a trapped atom 
interferometer in the simplest two-mode model. This will allow us to illustrate the basic 
physics. We then investigate how this simple two-mode model has to be modified when 
taking into account the many-body structure of the wave function. Optimal control 
techniques are applied to prepare the desired input states [311 l32l [33] . 

In our calculations we always assume zero temperature. The effects of additional 
dephasing and decoherence due to fundamental quantum noise and due to thermal 
excitations at finite temperature will be discussed in the last section. 

2. Two- mode model description of atom interferometry 

When atom interferometry is performed with a trapped Bose Einstein condensate 
(BEG), we consider the following key stages: The splitting stage (with the time duration 
^spiit), where the condensate wave function is split into two parts, the phase accumulation 
stage (Tphase)? where the atoms in one arm of the interferometer experience an interaction 
with some weak (classical) field, and finally the read-out stage (Ttof), where the phase 
accumulated is measured after the condensates have expanded in time-of-fiight (TOF). 

In our theoretical description, we start by introducing a simple but generic 
description scheme of an interferometer in terms of a two-mode model for the split 
condensate. Such a model has been also proven successful for the description of 
interference with spin squeezed states [34l[35] and of condensates in double wells p!8l[23] . 
To properly account for the many-boson wave function, we introduce the field operator 
in second-quantized form [36j 

i!{x) = ClL (t)L{x) + CiR (t)R{x) . (1) 

Here, a^j^ are bosonic field operators that create an atom in the left or right well, 
with wavefunctions 0L,i?(^), respectively. In many cases, we can properly describe 
the dynamics of the interacting many-boson system by means of a generic two-mode 
hamiltonian [f] [371 EHt |33] 

H = (^a[dj^ + 4a^) + K (a[a[a^a^ + a^^a^Raj^aj^^ . (2) 

Here Q describes a tunnelling process, that allows the atoms to hop between the two 
wells and k is the non-linear atom-atom interaction, that energetically penalizes states 
with a high atom-number imbalance between the left and right well. We treat as a 
free parameter, but choose a fixed ^ Uo/2. Uq is the eflFective one-dimensioanl (ID) 

i From now on we use conveniently scaled time, length, and energy units [33^, unless stated differently. 
First we set h = 1 and scale the energy and time according to a harmonic oscillator with confinement 
length aho = \/^/(j^^ho) — 1 and energy huji^o. Our considerations deal with Rb atoms, we 
therefore measure mass in units of the ^^Rb atom, and then time is measured in units of l/cj/^o = 1.37 
ms and energy in units of hbJ}^o = 27r • 116.26 Hz. 
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interaction strength along the direction where the potential is split into a double well. 
A more accurate way of how to relate n and [/q will be given in section |4} 

In the example we discuss in this manuscript the initial state of the interferometry 
sequence is prepared by deforming a confinement potential from a single to a double 
well. This corresponds to a situation where one starts with the two-mode hamiltonian of 
equation ^ for a large tunnel coupling, and then turns ofl[ fi, as a consequence 

of the reduced spatial overlap of the wavefunctions (/>L,i?(^)- 

In the beginning of the splitting sequence tunnelling dominates over the non-linear 
interaction, and all atoms reside in the bonding orbital (j)L{x) + (I)r{x). This results in 
a binomial atom number distribution. When Q is turned off sufficiently fast, and the 
dynamics due to the non-linear coupling plays no significant role, the orbitals (/>L,i?(^) 
become spatially separated, but the atom number distribution remains binomial. In 
contrast, when Q is turned off sufl&ciently slowly, the system can adiabatically follow 
the groundstate of the Hamiltonian ^ and ends up approximately in a Fock state. As 
we will discuss below, such states with reduced atom number ffuctuations are appealing 
for the purpose of atom interferometry. 

2.1. Pseudospin operators and Block sphere 

A convenient representation of the two-mode model for a many-boson system is in terms 
of angular-momentum operators [371 EHl EH] . Quite generally, the internal state of an 
ensemble of atoms which are allowed to occupy two states (here left and right) can be 
described as a collective (pseudo)spin J = Yl^=iji-> which is the sum of the individual 
spins of all atoms. Here the total angular momentum is and the projection m on 
the z-axis corresponds to states where, starting from a state where the left and right 
well are each populated with N/2 atoms, m atoms are promoted from the right to the 
left well. Through the Schwinger boson representation 

4 = ^ + ^r^l) ^ ^2/ = - ^r^l) ^ i^^^L^L - ^r^r) (3) 

we can establish a link between the field operators cll^r and the pseudo-spin operators. 
Jx promotes an atom from the left to the right well, or vice versa, and Jz gives half the 
atom number difference between the two wells. 

One can map the two-mode wave function onto the Bloch sphere [40], which provides 
an extremely useful visualization tool for the purpose of atom interferometry (see 
figure [T]). A state on the north pole corresponds to all atoms residing in the left well and 
a state on the south pole to all atoms in the right well. All atoms in the bonding orbital 
corresponds to a state localized around x = N/2. This is a product state where the atoms 
are totally uncorrelated. In this state the quantum noise is evenly distributed among 
AJy = AJz = ^/N /2, i.e., it has equal uncertainty in number difference (measured along 
the z-axis) and in the conjugate phase observable (measured around the equator of the 
sphere). Similar to optics with photons, or as discussed in the context of spin squeezing, 
quantum correlations can reduce the variance of one spin quadrature, for a given angle 0, 



CONTENTS 



6 




Figure 1. Bloch sphere representation of a number squeezed state with squeezing 
factor (defined in text) ^at ^ 0.2. AJ^ corresponds to number squeezing, and AJ^ is 
proportional to phase squeezing. In the ring below the sphere we show the polarization 
of the state along the x-axis, which is proportional to the coherence factor a = 2{Jx)/N. 
For squeezed states there is also a noise AJ^ in the polarization of the state on the 
equator. 



J(P = COS (/)Jz+^in (/)Jy at the cost of increasing the variance of the orthogonal quadrature: 
at the angle the variance AJ| becomes minimal, whereas the orthogonal variance 
AJ|^^^2 becomes maximal j41j- For example, the squeezed state shown in figure Tjhas 
reduced number fiuctuations, as described by the normalized number squeezing factor 
= ^Jz/ {^/N /2)^ and enhanced phase fiuctuations, as described by the normalized 
phase squeezing factor (^phase = ^Jy/{VN /2). 

Within the pseudo-spin framework, the two-mode Hamiltonian becomes 



H 



(4) 



which is completely analogous to the Josephson Hamiltonian of superconductivity. Q 
is associated with the (time-dependent) Josephson energy and the non-linearity tv with 
the charging energy [42j. For a given state on the Bloch sphere, the tunnel coupling 
rotates the state around the x-axis, whereas the non-linear part distorts the state such 
that the components above and below the equator become twisted to the right- and 
left-hand side, respectively. The twist rate due to increases with distance from the 
equator (see figure [2]) . 




Figure 2. Time evolution of states on Bloch sphere. The different panels report 
results for a (a) binomial, (b) number squeezed [^at ^ 0.2], and (c) phase squeezed 
state [^0 ~ 0.2]. The time interval is T = 12.8 for the phase squeezed state, and 
T = 32 otherwise. Due to the non-linear atom-atom interactions the states become 
distorted, thus spoiling the interferometer performance. 



2.2. Readout noise in the interference pattern 

We next address the question which states would be the best for the purpose of 
reading out an atom interferometer. For this purpose we neglect the non-linear atom- 
atom interaction that distorts the wave function during the phase accumulation time 
and postpone the question of how the specific states could be actually prepared in 
experiment. Our discussion (which closely follows Ref. [43j, with some extensions) is 
primarily intended to set the stage for the later discussion of the full atom interferometer 
sequence in presence of atom-atom interactions. 

We assume that in the phase accumulation stage the wave function has acquired 
a phase 9. Instead of describing the interaction process dynamically, which would 
correspond to a Bloch-sphere rotation of the state around the z-axis, we directly assign 
the phase to the single-particle wave function, such that the field operator reads [43j 

4f{x) = cll ^l[x) + ciR e-'^^R[x) . (5) 

To read out 0, one usually turns off the double well potential and lets the two atom 
clouds overlap in time of flight (TOF). The accumulated phase is then determined from 
the interference pattern [16j. After release, the wave function evolves under the free 
Hamiltonian. Here and in the rest of this paper we ignore the influence of atom-atom 
interactions during TOF which is justified in low dimensional systems [44j, but might 
be problematic under other circumstances [45l [46l [47j . 

As a representative example, we consider for the dispersing wavefunctions (/>l(^) 
and (j)R{x) two Gaussians with variance a^, which are initially separated by the 
interwell distance d. The density operator n{x) = ^'''(x)^^(x) associated with this TOF 
measurement can then be computed by using the pseudospin operators of equation ^ 



CONTENTS 



8 



as 

N / - -\ 

H^) = y (^^ + + ^V^lTIr [co8{kx + 9)Jx - 8m{kx + 9)Jyj + (n^ - nR)Jz . (6) 

In the foUowing we analyze the mean (ensemble-averaged) atomic density and 
its actuations around the mean value, such that for a single-shot mea- surement the 
outcome is with a high probability within the error bonds dened by the variance if the 
probability distribution is Gaussian. If the initial state preparation was for a symmetric 
double well potential, then the mean atom number difference vanishes ( J^) = 0, and for 
the real- valued tunnel coupling we can set ( J^) = 0. The mean density thus becomes 

n{x) = = y (^nL{x) + nR{x)^ + 2^s/nL{x)nR{x) cos{kx + 9) (^jj^ , (7) 

which is shown for representative examples in figure [3) Here, the visibility of the 
interference fringes [equation ([7|] is determined by the coherence factor 



a = 2{l)/N. (8) 



It is one for a coherent state with perfect polarization, (Jx) = where all atoms 

reside in the bonding orbital + e~^^(/)i^(x), while {Jy) = {Jz) = 0. In contrast, 

a Fock-type state, where half of the atoms reside in the left well and the other half 
in the right well, has no defined phase relation between the orbitals: it is completely 
delocalized around the equator of the Bloch sphere, and the coherence a = vanishes. 
In general, squeezed states have < a < 1. 

From equation ^ we can also obtain the density fiuctuations 

An(x)2 = {{n{x) - n{x)f) (9) 

= ^riLTiR ((a^') cos^kx + 9) + l^jf^ sm\kx + 9)^ + K " ^rY (^JI) . 

Noise contributions from all pseudospin operators contribute, diflFerently weighted by 
the time and space dependent probability distributions nL^R{x^t). The x-contribution, 
proportional to (AJ^), accounts for an uncertainty in the polarization of the state 
on the equator (see figure [T]). It is zero only for a coherent state. We will discuss 
the physical meaning of this quantity later in context of MCTDHB, section [4} The y- 
contribution accounts for the intrinsic phase width of the quantum state, which provides 
a fundamental limit for the phase measurement, and the z-contribution for the number 
fiuctuations between the two wells. 



2.3. Phase sensitivity 



The ideal state for detecting small variations is one where n{x) as a function of 9 
has a sufiiciently large derivative, and the fiuctuations An(x) are sufiiciently small. In 
order to resolve A0 the inequality 

dn{x) 



An{x) < 



09 



A9 



(10) 
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Squeezing factor 


Definition 


Useful squeezing 
Phase squeezing 
Number squeezing 


in = ^Jy/{a^/Nl2) 
= AJ,/(^/2) 



Table 1. Definition of squeezing factors used in this work. is the useful squeezing 
that determines the sensitivity of an interferometer. It differs from the phase squeezing 
through the coherence factor a = 2{Jx)/N. is the number squeezing that depends 
on the number fluctuations between the two wells. 



has to be fulfilled. The explicit expression in terms of Ul^ji follows immediately from 
equations ([T]) and ([9]). A particularly simple expression is obtained if we keep in 



equation ( 10) only the dominant contribution from the phase noise, which is proportional 
to {Jy) (other contributions depend on nL^R{x) and are expected to be less important), 
which leads to the estimate for the phase sensitivity in terms of shot noise by Kitagawa 
and Ueda [H] (see also table [T]) 



iR=^ -. (Ill 



(4) 

This expression will serve us as a guiding principle for optimizing atom interferometry. 
We will consider in the following the phase squeezing at the optimal working point. 
The normalization by {Jx) takes into account that improving interferometric sensitivity 
requires not only to reduce noise but also to maintain a high interferometer contrast a 



(which determines the diflFerence to ^phase via = Cphase/<^) • Equation (11) provides a 
simple way to estimate the sensitivity obtainable by a given initial state. For a coherent 
state with a binomial atom number distribution, the mean value and the variance are 
proportional to the total atom number A^. Thus, ^r = 1 and the phase sensitivity 
A9 = is shot-noise limited. We refer to this limit as the standard quantum limit. 

For the purpose of reading the interference pattern, phase squeezed states are ideal 
because they allow to reduce the sensitivity below shot noise, ^r provides a measure 
of useful squeezing for metrology ^48j: a state with < 1 allows to overcome the 
standard quantum limit by a factor ^r. The lower bound of the sensitivity is provided 
by the Heisenberg limit. From the commutation relation [J^, Jz] = iJx one obtains 
the uncertainty relation AJyAJ^ > N/2. Thus, for a state with a maximal number 
uncertainty, AJ^ = A/2, the standard deviation AJy is on the order of unity, and the 
Heisenberg limit becomes ^r = ^2/A. 

In figure [3] we show density and density fiuctuations for a (a) binomial, (b) phase- 
squeezed, and (c) number squeezed state. One observes that the fiuctuations are smallest 
for the phase-squeezed state, and become larger for the coherent and number squeezed 
state. To appreciate the sensitivity of the interferometer, we also plot the density profile 
for a state that has acquired a phase 9 of the order of shot noise. As apparent from 
the insets, which magnify the regions of smallest noise, the smallest 9 variations can be 




Figure 3. Density and density noise for (a) binomial, (b) phase squeezed, and (c) 
number squeezed state, and for a phase = (green areas) and = 1/y/N (red 
areas), corresponding to shot noise. We use N = 100, width a = 0.05, a read-out time 
Ttof = 10 and an interweh separation ofd = 5. The degree of squeezing is characterized 
by a squeezing factor of ^at = 0.18 (^^ = 0.18). The dashed hues show the density 
profile of the condensates without interference. The insets magnify the region where 
the noise is least. The phase sensitivity is best for the phase squeezed state, although 
in certain regions noise is enhanced due to AJx. 



resolved with the phase-squeezed state. 

In a typical double well interferometer experiment, where the interference is read 
out in TOF, the exact number of atoms is not known, and one has to obtain and 
the accumulated phase 9 from a suitable fitting procedure. To benefit from the regions 
of reduced density noise, the density distribution n{x) has to be weighted appropriately 
by An{x). 

The dependence on the duration of the read-out stage is shown in figure [4} For short 
Ttof the overlap between the condensates is small, and the contribution from AJ^ (the 
third term in equation ([9])) important. Since small phase fluctuations are accompanied 
by large number fluctuations, there exists an optimal degree of phase squeezing for the 
interferometer input state. For large enough Ttof this effect is less important and 
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Figure 4. (a) Best achievable versus read-out time Ttof foi" atom number N = 100, 
a = 0.05, and d = 5 (as in Fig. [s]). For given Ttof, the optimal phase squeezing ^phase 
(dashed line) is found such that (solid black line) is minimized. approaches the 
Heisenberg limit (gray line) for large Ttof. The difference between and ^phase is due 
to a < 1 (b) Polarization noise, AJx, of the corresponding phase squeezed states of 
(a). 
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approaches the Heisenberg hmit. We optimize the degree of phase squeezing of the 
initial state (i.e., at Ttof = 0) for given Ttof, such that is minimized. This is shown 
in (a) versus Ttof (sohd hue), together with the corresponding optimal phase squeezing 
<^phase (dashed line). For small ^phase, polarization noise AJ^^ [shown in panel (b)] is 
introduced. This noise reduces the spatial region where the sensitivity is high [see also 
figure [3](b)]. In conclusion, the expansion in TOF has to be sufficiently long, such that 
the atom clouds can expand sufficiently far, and the phase sensitivity is no longer limited 
by number ffuct nations between the clouds. 



2.4' Interferometry in presence of atom- atom interactions during the phase 
accumulation stage 

In the following we discuss how the atom-atom interactions affect the phase 
distribution of the split condensate during the phase accumulation stage and spoil atom 
interferometry, and what could be done to minimize those effects. 

For the purpose of atom interferometry, it is convenient to introduce the phase 
eigenstates [19] 

^^e^*-|m), (12) 



2tt 

Hi 

where |m) is a state with atom number imbalance m between the left and right well. 
The relation between \m) and \(f)) corresponds to a Fourier transformation. We can now 
project any state on the phase eigenstates and obtain the phase representation of the 
state. From this we obtain the phase width A(/) (see also [49j). The noise AJ^/, which 



enters the phase sensitivity of the TOF-interferometer according to equation (11), has 
a similar time evolution as the phase width . However, while the phase width is bound 
to values below 27r, the upper bound of A is given by the total atom number A^, and 
therefore AJ^ = yfN/\(j). 

Phase diffusion due to atom-atom interactions is best illustrated in the Bloch sphere 
representation (see figure [2]): The non-linear coupling kJ^ from equation ^ twists the 
state, and the larger the number fiuctuations AJ^ the faster the state winds around 
the Bloch sphere. In the notion of phase eigenstates, a state with a well defined phase 
has a broad atom- number distribution. As each atom- number eigenstate \m) evolves 
with a different frequency, the time evolution of a superposition state will suffer phase 
diffusion. The phase width broadens with rate \T[\ 

R = SUoAn = 8UoCn{Vn/2) . (13) 

As a result, states with small number ffuctuations {^n < 1) are more stable during the 
phase accumulation time. 

One immediately sees a conffict of requirements. For best readout of the interference 
pattern we want phase squeezed states, but those are very fragile and result in a fast 
phase diffusion and a short measurement time. On the other hand, number-squeezed 
states allow for longer measurement times but have a rather poor readout performance. 
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Figure 5. Example for the phase sensitivity of a binomial (solid line), a number 
squeezed (dashed line), and a phase squeezed (dashed-dotted line) state with squeezing 
factors ^ 0.22 and ~ 0.22, respectively. Parameters are N = 100 and U^N = 1. 
The binomial and phase squeezed states have a better initial phase sensitivity, whereas 
the number squeezed is much more stable against phase diffusion, and has a much 
better sensitivity at later times. The gray line corresponds to shot noise {^r = 1). 



In the remaining part of the paper we wiU discuss the optimal strategy for interference 
experiments, and how one can implement them in realistic settings. 

To demonstrate the above reasoning, we show in figure [5] the achievable phase 
sensitivity as a function of phase accumulation time Tphase in presence of atom-atom 
interactions for a binomial (solid), a number squeezed (dashed line), and a phase 
squeezed (dashed-dotted line) state with squeezing factors ^ 0.22 and ^ 0.22, 
respectively. The phase squeezed state has initially sub-shot noise phase sensitivity. 
For longer hold times, the number squeezed state outperforms the phase squeezed and 
binomial ones due to its smaller phase diffusion rate. 



3. Optimizing atom interferometry 

When designing a trapped atom interferometer for measurements, one has to consider 
the conflicting requirements from phase diffusion and readout, the one asking for number 
squeezing, the other for phase squeezing. In the following we will outline a few strategies 
of how to optimize interferometer performance for realistic double well settings. 



3.1. A very simple estimate 

We first analyze the impacts of trap geometry and initial state preparation. We start by 
employing a simple model to illustrate the effects of the non-linearity in the evolution 
of the split trapped BEC. Let us first look at the interaction energy of a trapped cloud, 
and how it changes with the number of trapped particles. Without loosing generality, 
we discuss in the context of a harmonic trap characterized by the mean confinement ujq 
and the length scale ao, which are defined through 

ujo = ^uj^^UyUJ^ , ao = y — ^ . (14) 
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With N atoms in the trap the chemical potential in Thomas-Fermi approximation 
[HQ] is given by 

^ = — \-^) ' ^''^ 
where is the 5- wave scattering length. Adding a single atom to the trap, changes ji 
by 

dN = lN- ^'^^ 
This quantity corresponds to the effective ID interaction parameter [/q discussed earlier. 
We can now estimate the scaling of the phase diffusion rate R caused by a number 
distribution with fluctuations AA^ = ^nVN /2 after splitting = 1 corresponds to a 
binomial number distribution, whereas < 1 to number squeezing). 



2 



R(x^N A"To (jj^ al ms . (17) 

1 ^ 

For ^''Rb atoms (a^ = 5.2 nm, m = 87) one obtains R ^ 0.022 A~ To cJq s~-^, or in 
scaled units R ^ 0.29 N lo cjq . If we now set the phase diffusion rate R equal to 
the phase accumulation rate (the signal we want to measure), we find the limit 

^AE^^^ = R for the sensitivity of a single-shot interferometer measurement, even for 
perfect readout. It is interesting to note that this sensitivity limit is only very weakly 
dependent on the atom number N. 

If the interferometer measurement is not limited by readout, we can identify the 
following strategies for improving the interferometer performance (or, equivalent ly, 
reducing the effect of phase diffusion). 

Minimize the scattering length. The best is to set = 0, which can in principle 
be achieved by employing Feshbach resonances [51j. Drastic reduction of 
phase diffusion when bringing the scattering length close to zero was recently 
demonstrated in two experiments [52l [5^ . A disadvantage thereby is that using 
Feshbach resonances requires specific atoms and specific atomic states. These states 
need to be tunable, and are therefore not the 'clock' states usually used in precision 
experiments which are immune to external disturbances like magnetic fields. 

Choose a trap with weak confinement. This route seems problematic, since the 
timescale in splitting and manipulating the trapped atoms scales with the trap 
confinement. Optimal control techniques, like discussed in Refs. jSU |32[ |33] , will 
be needed to allow splitting much faster than the phase diffusion time scale. It 
is interesting to note that one needs a strong confinement only in the splitting 
direction. In the other two space directions the confinement can be considerably 
weaker. This suggests to work with strongly anisotropic traps. 
For an elongated cigar-shaped trap (ID geometry) with confinement ratio Cid = ^ 
(strong confinement (jJ± in radial directions, weak confinement in axial direction) 
one finds 

i'""^ = UJ^VC^ ■ (18) 
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For a flat pan cake-shaped trap (2D geometry) with confinement ratio C2D = '^^^^^ 
(strong transverse confinement (jJ± and weak in-plane confinement cjpiane) one finds 




With a confinement ratio C ^ 1/1000, which is easily obtainable in experiments 
the phase diflFusion is reduced by a factor 10 in a ID geometry and a factor 100 in 
a 2D geometry. 

Increase number squeezing in the splitting process. This directly reduces the 
phase diflFusion rate and hence leads to a better limit for the minimal detectable 
signal. Number squeezing can be achieved during the splitting process. It is 
mediated by the atom interactions and one has to achieve a careful balance between 
the interactions necessary to obtain sizable number squeezing and the decremental 
effect of the interactions during the phase accumulation time. This will be one of 
the central parts in our optimization discussed below. 

In an ideal interferometer one would like to use clock states, create strong squeezing 
during the splitting process, exploiting the non-linearity in the time evolution, and then 
turn off the interactions (by setting the scattering length to = 0) after splitting. All 
together might, however, be difficult or even impossible to achieve. In the remainder of 
the manuscript we will discuss the different contributions to the precision of an atom 
interferometer, and investigate how the performance can be optimized. 



3.2. Optimization of the many-boson states 

3.2.1. Optimized Number Squeezing First we shortly discuss how number-squeezed 
states can be prepared during the splitting stage. Number squeezed states are created 
in a double well potential when the interaction energy starts to dominate over the 
tunnel coupling, the latter being controlled by the barrier height and the double well 
separation. A natural way to achieve high number squeezing is dynamic splitting of 
a BEC |49l [54j, such that the wavefunction can adiabatically follow the ground state 
[38j. However, this may take very long, possibly longer than the phase diffusion time 
of the split condensate. In our earlier work [|32[ [33l |39] we employed optimal control 
theory (OCT) to find splitting protocols which allow for high number squeezing on 
a fast timescale, at least one order of magnitude shorter than for the quasi-adiabatic 
splitting. These protocols can be viewed as the continuous transformation of a (close to) 
harmonic potential into a double well. Thereby, a barrier is ramped up at the center, and 
simultaneously the two emerging wells are separated. In many cases this splitting process 
can be parameterized by a single parameter, whose time variation is obtained within 
optimal control theory such that sizeable squeezing is created, condensate oscillations 
are prevented after splitting, and phase coherence is better preserved at the end of the 
splitting process [32l [33l [39j . 

Unless stated otherwise, in the following discussion of the dynamics during the 
phase accumulation stage we use number-squeezed states as initial states which are 



CONTENTS 



15 









20 40 60 

Phase accumulation time 



80 



Figure 6. (a) Optimal phase sensitivity versus phase accumulation time for N = 100 
(solid lines), N = 500 (dashed lines), N = 2000 (dashed-dotted lines), and TV = 8000 
(dots). Interaction strength is such that UqN = 1. The black lines show results for 
number-squeezed states, and the gray lines result for slightly tilted number-squeezed 
states, (b) The corresponding optimal number squeezing, (c) The total tilt angle 
0tiit = ^tiit^puise is determined by the pulse duration (here Tp^ise = 2) and strength 
^tiit (very similar for all N). The diamond symbols in (a) and (b) show estimates from 
a simple model. 



obtained as the ground states of equation (|4]) for finite values of tunneling Q. They are 
very similar to those obtained by OCT the splitting. Similar initial states obtained by 
exponential splitting have a smaller degree of coherence. During the phase accumulation 
stage, we set = 0. 

How much number squeezing is ideal for a pre-determined phase accumulation time 
of an interferometer? In figure [6] we show results where we optimize the degree of number 
squeezing for the interferometer input states at Tphase = 0, in order to achieve the best 
phase sensitivity at a given time Tphase- The top panel reports the best achievable phase 
sensitivity for an initially number squeezed state (black lines), and the middle panel 
reports the corresponding number squeezing. For short phase accumulation times, less 
initial number squeezing is better. With increasing accumulation time more number 
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squeezing becomes favorable. This is due to the competition between phase fluctuations 
AJ^, which increase with number squeezing, and the decrease of phase diffusion for 
states with high number squeezing. 

The results can be well explained by a simple model. Neglecting the effects of 
reduced phase coherence, we have initially = Cphase' initial phase squeezing. 



Phase diffusion with rate R then results in = y (Cphase)^ + ^^^phase- i^ey^t use 
that ^A^Cphase ^ 1' which is in the spirit of the Heisenberg uncertainty principle and 
agrees well with our OCT results. Putting in all the constants, we have 



in = ]^ ^ + ^m%U^T^^^^. (20) 

The minimum with respect to is found as = 1/ {2^yU^NT~y^)^ which yields a 
best phase sensitivity ^™ = 2^2[/oA^Tphase- For a given flnal we see that Tphase is 
indirectly proportional to the interaction parameter [/q, which allows rescaling of Tphase 
in case of a different Uq. 

Predictions of the simple model are shown in flgure [6] by the diamond symbols. 
The agreement with the exact results is very good in (a), and gives the right scaling 
in (b). Indeed, ^™ in flgure [g] is independent of for long times, as long as UqN is 
constant. This is not true for <^]^^^. The reason is that neglecting phase coherence makes 
the minima with respect to much more shallow. For longer times, however, phase 
coherence becomes more important and the present approximations are no longer valid 
for small N. 



3.2.2. Optimized specialized initial states A different strategy for improving the 
interferometry performance is to prepare the system at the beginning of the phase 
accumulation stage (i.e., at Tphase = 0) in a special state, which evolves under the 
influence of the non-linear interaction after some pre-determined time into a state with 
high intrinsic sensitivity. The ideal initial state would be the time reversal of a phase 
squeezed state. We denote this strategy as refocusing. When the condensate is released 
at the optimal time, and expands in absence of interactions to form the interference 
pattern, interferometry can be performed with a sensitivity determined by the properties 
of the refocused state. This can be achieved because the phase accumulation (rotation 
around z-axis on the Bloch sphere) and the non-linear coupling do not interfere. 

Such a refocusing strategy is related to spin echo techniques, which where 
investigated by turning the scattering length from repulsive to attractive [55j. 
However, the latter has given a rather poor improvement, because it does not lead 
to a perfect time reversal of the many-body dynamics [56j . Artiflcial preparation of the 
desired time reversed states, seems to be very diflicult. We did not succeed in this task 
using optimal control techniques. 

One state that leads to very good refocusing can be prepared by tilting the initial 
number squeezed state on the Bloch sphere slightly against the direction of the twist 
originating from J^. The tilt can be achieved by applying a short tunnel pulse within a 
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Figure 7. Refocusing control sequence for N = 100, UqN = 1, Tp^ise = 2, and 
(a) Tphase = 1 as well as (b) Tphase = 4. The lower panel shows ^i^, while the upper 
panel the coherence factor a. In the first stage, the tunnel pulse is applied and tilting 
is achieved. In the second stage, the state refocuses to a state with a good phase 
sensitivity, in (a) below shot-noise. The Bloch spheres visualize the time evolution. 

time interval Tp^ise that rotates the number squeezed state. In real space, this operation 
corresponds to lowering the barrier for a suitable amount of time, which cannot be 
done arbitrarily fast because of condensate oscillations. Appropriate controls of the 
barrier will be discussed in detail in context of a realistic modeling in section [4} Within 
the generic model we consider for simplicity square-f] pulses, which is the best possible 
pulse in presence of interactions [57j. Examples are shown in figure [T] for different Tphase- 
During the tilting pulse sequence and the phase accumulation stage, 'rephasing' happens 
and the phase fluctuations decrease. Simultaneously, the phase coherence is restored to 
a value close to one. This significantly improves the phase sensitivity, for short times one 
can even reach below shot noise. However, this cannot be done perfectly, the degree of 
phase squeezing <^phase achieved after refocusing is always less than the degree of number 
squeezing of the original state. 

We next optimize systematically both parameters of initial number squeezing and 
tilt angle. The lowest panel of figure [6] shows the optimal tilt angles. From the upper 
panels (bright lines) we find a clear improvement of phase sensitivity for a given phase 
accumulation time. The dependence on is very distinct now for small atom numbers, 
and saturates for large N. We find that for small N the improvement is roughly a factor 
of three in time, and for large N approximately an order of magnitude. 

3.3. Optimize trapping potential and atom number, esults of generic two-mode model 

We next proceed to a more detailed analysis of the ideal trap parameters. Considering 
our previous discussion of section |3.1[ we expect an improvement of the interferometer 
performance when increasing the anisotropy of the trap. In the following calculations 
we choose a fixed N = 100 (the results are not expected to depend decisively on N). 
We first investigate the role of the interaction UqN^ and the pulse duration Tp^ise for 
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Figure 8. Best achievable phase sensitivity versus interaction strength for various 
pulse durations for (a) Tphase = 1 and (b) Tphase = 20. Atom number is = 100. We 
optimize for the initial number squeezing and the tunnel pulse. 



the refocusing strategy. In figure [8] we plot the best phase sensitivity using refocusing 
versus interaction strength for (a) Tphase = 1 and (b) Tphase = 20. Sub-shot noise 
phase sensitivity is clearly achievable for short Tp^ise, or for U{)N <C 1. Short pulses 
are favorable and give better phase sensitivity. This is in particular important for short 
phase accumulation times Tphase (see figure |6]) . Optimizing the pulse form and duration 
within a realistic modeling of tilting on the Bloch sphere will be discussed in detail in 
section HI 

Quite generally, we can expect that for a reduced interaction parameter it is more 
difl&cult to obtain high number squeezing in the splitting process. To estimate the time 
scale for achieving a certain degree of number squeezing, we consider splitting protocols 
derived in a previous work within the framework of optimal control theory [33], and 
discussed already in section [3.2. 1[ They directly provide us optimized number squeezing 
for a given splitting time Tgpht. For a given phase accumulation time Tphase, we optimize 
then the tunnel pulse which tilts the number squeezed state, similar to the analysis of 
section 3.2.2, The best possible phase sensitivity for various Tgpht and U{)N values is 



shown in figure [9] for (a) Tphase = 1 and (b) Tphase = 20. For both cases sensitivity 
distinctly sub-shot noise sensitivity can be achieved. To achieve the squeezing needed 
to boost interferometer performance, a finite [/q is needed, and for a given Tgpht there 
exists an optimal value of [/q. This value decreases for longer splitting times. 

In order to analyze the dependence on A^, we consider a realistic 3D cigar-shaped 
trap with transverse trapping frequency cj^ = 27r x 2 kHz as typically realized in atom 
chip interference experiments. The effective ID interaction strength in the splitting 
direction [/q is then approximately proportional to C^^^ . A more rigorous estimate, that 
is used in the calculations, is given in Ref. [33]. For = 100, U{)N ^ 1(0.1) corresponds 
to an aspect ratio of Cid ^ 1/100(1/1000). These values change for A^ = 1000 to 
CiD ^ 1/1000(1/10000). For the pan cake-shaped trap we have C2D ^ a/Cld, and 
thus U^N = 0.01 is within reach for C2D ^ 1/1000 and A^ = 100, C2D ^ 1/10000 and 
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Figure 9. (a,b) Phase sensitivity for a sequence including splitting and phase 
accumulation, versus interaction strength U{)N and Tsput- Parameters are N = 100, 
(a) Tphase = 1, (b) Tphase = 20, and Tp^ise = 2. For a given Tsput we take the best 
number squeezing achieved by OCT. Number squeezing is higher for larger Tgpiit and 
Uq values. We also optimize for a tunnel pulse which rotates the number squeezed 
state. The minimum phase sensitivity decreases very slowly with Tgpiit . The contour 
lines are at (a) = 0.35 and (b) = 0.6. 



= 1000. 

Let us first consider the case without refocusing. We find approximately ^ 
2y^2[7oiV7^hase, and, considering the dependence of Uq on the trapping potential, we 
obtain 

- al/'uTcl/^N'/Xil . (21) 

In order to reach sub-shot noise phase sensitivity we need the confinement ratio 
CiD, atom number A^, and Tphase all to be small. In figure 10 (a) ^r is plotted for 
CiD = 1/100, 1/300 and 1/1000 (solid, dashed, and dashed-dotted line, respectively). 
As we have seen in section [3^ equation (21) is valid only if the coherence is well 

^2 

preserved. We estimate breakdown of this approximation when a ^ 1 — ^ ^ 0.6. 
From this we can obtain the time after which ^r is expected to grow rapidly because 
the coherence factor tends to zero. It is given as TJ-oh — , 2/5 6/5^6/15 5 and shown in 

2" 15 ' CLg ^ J_ ID 

figure 10 (f) for different Cid- 

We now turn to refocusing. In the ID elongated trapping geometry we find a more 



moderate increase of ^r with N compared to the case without refocusing, see figures 10 



(a) and (b). This illustrates that the refocusing works better for large N, as long as 
UqN is constant, see also figure |6) 

The accumulated phase for a potential AE is given as ^ = AE'Tphase, and the 
smallest detectable potential difference in time Tphase becomes AE""^^^ = ^r/ (VNT^hase)- 
Without refocusing, we find an improvement with N and Tphase, AE"™ ~ ^~^^^'^phlse^ 
figure 10 (c), (d) and (e). Similar scalings also holds for the case with refocusing. We 
expect that for 2D traps, AE""^^^ somewhat below 10~^ is within reach. 
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Figure 10. Optimization for a realistic cigar shaped trapping potential with transverse 
frequency uo^ = 27r • 2 kHz (cj^ ^ 17 in scaled units) and aspect ratio Cid = 1/100 
(solid lines), Cid = 1/300 (dashed lines), and Cid = 1/1000 (dashed-dotted lines), 
(a), (b) versus TV for (a) Tphase = 1 and (b) Tphase = 20. (c), (d) Minimal detectable 
potential AE"^^" versus TV for (c) Tphase = 1 and (d) Tphase = 20. (e) AT;"^^" versus 
time for = 1000 (the vertical dashed lines in (a)-(d) denote cuts at = 1000). 
The black lines show results for the case when only the initial number squeezing is 
optimized. The coloured lines show results with refocusing, i.e., where both initial 
number squeezing and tunnel pulse are optimized, (f) Time Tcoh after which is 
expected to grow rapidly due to loss of coherence, if refocusing is not applied. 



This confirms, in agreement with the scaling analysis of section [3^11 that a stronger 
trap anisotropy appreciably helps to reduce phase diffusion and, in turn, to improve 
the phase sensitivity of the interferometer. Possible limitations of strongly anisotropic 
systems are discussed in Sec. [5) The atom number helps to improve absolute 
sensitivity, but makes it more difficult to demonstrate measurements with a sensitivity 
below shot noise. 

4. Interferometer performance within MCTDHB 

Until now we have used a generic two-mode model to describe the interferometer, which 
captures the basic processes and physics, but ignores many details of the condensate 
dynamics in realistic microtraps. More specifically, the modeling of the splitting 
process and of the rotation pulses requires in many cases a more complete dynamical 
description in terms of the multi-configurational time dependent Hartree for Bosons 
method (MCTDHB) [58j. In this section we discuss first the MCTDHB details relevant 
for our analysis, and its relation to the generic model. The main part will be concerned 
with the simulation and optimization of tunnel pulses for achieving tilted squeezed 



states, as discussed in section 3.2.2 in context of refocusing. An exhaustive discussion 



of MCTDHB, as well as optimal condensate splitting can be found elsewhere 
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In the two mode Hamiltonian of equation ([2]) we did not explicitly consider the 
shape of the two orbitals (j)L and (j)R^ but lumped them into the effective parameters 
Vt and K.. The dynamics is then completely governed by the wavefunction accounting 
for the atom number dynamics. Within MCTDHB, both the orbitals and the number 
distribution are determined self-consistent ly from a set of coupled equations, which are 
obtained from a variational principle. This leads us to a more complete description, 
accounting for the full condensates' motion in the trap. The state of the system is 
then given by a superposition of symmetrized states (permanents), which comprise the 
time dependent orbitals. Instead of left and right orbitals, such as used in the two- 
mode model, we employ for the symmetric confinement potential of our present concern 
orbitals with gerade and ungerade symmetry. The time dependent orbitals then obey 
non-linear equations, which depend on the one- and two-particle reduced densities [59j 
describing the mean value and variances of the number distributions [60j . 

The atom number part of the wavefunction obeys a Schrodinger equation with the 
Hamiltonian 

U = + ^ XI ^Wq^l^rnWkqlm , (22) 

k,q,l,m 

where the indices are either g (gerade) or e (ungerade or excited). We observe that, 
in contrast to the two-mode Hamiltonian of equation ([2]), the atom-atom interaction 
elements 

Wkqlm = Uo J dx(f)l(x,t)(f)l(x,t)(l)i(x,t)(f)m(x,t) , (23) 

as well as the tunnel coupling Q = J dx(l)l(x)h(j)e(x) — J dx(j)''g(x)h(j)g(x) are governed 
by the orbitals. The only input parameter of the MCTDHB approach is the trapping 
potential V\(x)^ which enters in the single-particle Hamiltonian h(x) = — (V^/2)+Va(^)- 
We note that Q obtained within MCTDHB cannot be directly interpreted as tunnel 
rate, but has to be renormalized if the two-body matrix elements diflFer from each 
other [38l [6T]. Thus, there is in general no direct correspondence between the two- 
mode model and the MCTDHB approach. MCTDHB, which relies on time-dependent 
orbitals, captures a large class of excitations not included in a two-mode model. If the 
calculations converge when using more modes, MCTDHB reproduces the exact quantum 
dynamics, as discussed in [62], ?. I63]. 

In our MCTDHB calculations we consider a cigar-shaped magnetic confinement 
potential prototypical for atom chips ^ . Splitting is assumed to be along a transverse 
direction and is accomplished using rf dressing [GU [65l |66l [G^ • For illustration purpose, 
the trapping potential in splitting direction can to very high accuracy be described by 
a quartic potential of the form: 

V(x,t) = a(t)x^ + b(t)x^ , (24) 

where for most cases b(t) varies very slowly and can be assumed as constant. This 
potential grasps the essential features of the initial and final potential and the time 
evolution a(t) describes how the potential is split and the barrier is ramped up. a(t) 
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large and positive characterize the initial single well, a{t) large and negative the split 
double well, the constant h{t) the confinement during the splitting. In the calculations 
we use the exact form of the potential used in atom chip double well experiments [64j . 
To describe the transformation of the potential we introduce a control parameter \{t) 
connected to the amplitude (phase) of the RF field. Thereby, values of — 2/3 < A < 
translate to a single well, and A ^ 1 to a double well. 

Within MCTDHB the pseudospin operator Jx has to be rewritten in terms of 
the gerade and ungerade orbitals. In the new basis it measures the atom number 
difference with respect to the two states, Jx = — ^^^g), as discussed in more 

detail in the appendix. It is important to note that the gerade and ungerade orbitals 
are natural orbitals, i.e., they diagonalize the one-body reduced density of the system 
[59]. Therefore, if both of them are macroscopically populated one obtains a fragmented 
condensate [36j. We can thus interpret the coherence factor a as the degree of 
fragmentation. The system has maximal coherence, if its state is not fragmented 
but forms a single condensate. Coherence is lost if the condensate fragments into two 
independent condensates. In between, we have a finite, but reduced coherence. Similarly, 
we can interpret AJx as the number uncertainty between the fragmented parts. 

Optimal control theory (OCT) |68[|35, is a very powerful tool to find a path which 
optimizes for a certain control target. In our earlier work [32l|33] we implemented and 
optimized condensate splitting within MCTDHB [68l |3T]. We found that, although 
the generic two-mode model describes qualitatively the splitting dynamics, the more 
complete MCTDHB description is needed for a realistic modeling. This is because 
one needs to control condensate oscillations during splitting, and to ensure a proper 
decoupling of the condensates at the end of the control sequence. In this section we 
employ OCT to find the appropriate paths in varying the trapping potential to achieve 
the desired tilting on the Bloch sphere in short time Tp^ise and prevent excitation of 
condensate oscillations. 



Parameter correspondence between the models 

For the calculation of time-dependent condensate dynamics including oscillations, a self- 
consistent approach like MCTDHB is mandatory. However, we expect the generic model 
(equation ^) to properly describe the phase accumulation stage (f] = 0), provided the 
condensates are at rest. The two-body overlap integrals of the orbitals from MCHB 



(time-independent version of MCTDHB ^GOj), given in equations (23), are then constant 
and coincide. This is because (j)g and (/)e have degenerate moduli for split condensates. 
When comparing equations ^ and (22), we find the value of tv to be used in the generic 
model. Similarly, in the context of optimized splitting protocols in our earlier work 
[321 [33] . we have found that both optimizing Q{t) in the generic model and optimizing 
X{t) within MCTDHB yield the same amount of number squeezing for a given time 
interval. 



In figure 11 we show how the two-body overlap integrals from MCHB vary with 
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X(G) 



Figure 11. Two-body matrix elements for the ground states of a magnetic trap [64] 
versus splitting distance for N = 100 and UqN = 1. The splitting is parametrized 
by the parameter A. Wgege and Wggee coincide for the ground states, however not in 



general, see the examples of figures 12 and 13 



the control parameter that determines the shape of the confinement potential. During 
the transition from a single-well to a double well, they drop by roughly a factor of 
two. This is because in the final state the gerade and ungerade orbitals are delocalized 
over both wells. After reaching a minimum around A ^ 0.7, the overlap integrals start 
to increase again slightly. In context of splitting we found that it is reasonable to 
assume = [/o/2 throughout the splitting process [33], i.e., to take the value in the 
most relevant regime during condensate breakup (A ^ 0.7 — 0.8). Atom interferometry 
has to be performed with split condensates. This requires a splitting distance of some 
micrometers, corresponding to A > 1. With the corresponding value oi n 0.65[/o, 
phase diffusion in the phase accumulation stage can be well described using the generic 
model. This has in particular the advantage that we can translate our findings for the 



optimal states for atom interferometry from section |3.2.2| to MCTDHB calculations of 
tunnel pulses. 

4.2. Pulse optimization 

A key ingredient in interferometer performance is the preparation of an optimized initial 



state, as discussed in section [3. 2. 2[ This can be achieved by a 'tunnel pulses' facilitating 
the tilting of the initial number squeezed state. This operation has been analyzed by 
Pezze et al. in the context of a cold atom beam splitter. They used the generic model 
[57] and studied to which extent the creation of phase squeezed states from number- 
squeezed states is spoiled by atom- atom interactions. The real space dynamics has been 
neglected. 

In real space, the confining potential has to be modified to bring the condensates 
together for the tunnel pulse which accomplishes the desired 'tilt' on the Bloch sphere. 



As has been discussed in section 3.3, the duration of the tunnel pulse (Tp^ise) is very 
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critical for the interferometer performance, the pulse should be as short as possible. 
This has to be done without significantly disturbing the many-body wave function. To 
design appropriate control schemes with the shortest possible time duration we will have 
to take the real space dynamics of the BEC into account using MCTDHB. 

To find a control sequence for the preparation of the optimal initial states at 
^phase = for a givcu interferometer sequence within MCTDHB, we first start, following 



section 4.1, with the optimized initial number squeezing and the tilt angle required for 



rephasing at a given Tphase as calculated in section 3.2.2 within the generic model. We 



then choose the top of a Gaussian shaped initial guess 

\{t) = Ao + (Aend - Ao) • t/Tp^ise " A [g-^^^f^ - . (25) 



(^~^pulse/2)^ (^pulse/2)^ 



We start with a state of the double well potential with the required initial number 
squeezing. Then, as A decreases, the barrier is ramped down, the condensates approach 
each other. The desired AJ^ to be reached at Tp^ise is fixed by the required tilt angle. 
It can be tuned by the parameters A and S, corresponding to the depth and the width 
of the control parameter deformation, respectively. Finally, a double well V\^^^ is re- 
established, which completely suppresses tunneling of the two final condensates, at least 
if they are in the ground state. 



Results of our MCTDHB calculations are shown in figures [12] and [13] for interactions 
U^N = 0.1 and [/qA^ = 1, respectively. The pulse achieves the desired [dashed lines in 
(d)], as we expected from the generic model (dashed-dotted lines). However, it not only 
affects the atom number distribution, but also leads to an oscillation of the condensates 



in the microtrap. This can be seen in the density, which is depicted in figures 12 [^b) 



and 13(b). Condensate oscillations during the phase accumulation and release stage are 
expected to substantially degrade the interferometer performance, and may even lead 
to unwanted condensate excitations [63]. These oscillations can be avoided by using 
more refined tilting pulses, which can be obtained within the framework of OCT, where 
we now optimize for phase squeezing and a desired A at the final time of the control 
interval Tp^ise, corresponding to Tphase = 0. Some details of our approach are given in 
the appendix. 

For weak interactions a short tunnel pulse Tp^ise = 2 can be easily found which 
properly tilts the atom number distribution and brings the condensates to a stationary 
state at the end of the process. A typical control sequence for U{)N = 0.1 and optimized 
for a short Tphase = 1 is shown in figure [T2](a) [green solid line], the corresponding density 
given in (c). In panel (d) we depict the phase sensitivity (solid line), which compares 
in the phase accumulation stage very well with the desired behavior given by the generic 
model (dashed-dotted line). 

The gain of our OCT solution with respect to an 'adiabatic' control of Gaussian type 



[equation (25)], where A is modified sufl&ciently slowly in order to suppress condensate 
oscillations, depends on the chosen value of Tphase- Smaller Tphase values require larger 
tunnel pulses [see figure |6](c)]. In our example, the OCT pulses can be at least one 
order of magnitude faster, which means an improvement of up to 30% in ^r [compare 
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Figure 12. (a) Typical OCT-control (green solid line) for weak interactions UqN = 
0.1, atom number N = 100, Tp^ise = 2 and Tphase = 1^ compared to a simpler control 
(blue dashed line). The dashed vertical line separates the control sequence from the 
phase accumulation time thereafter. The corresponding condensates density are shown 
(b) for the initial guess and (c) for the OCT-solution, and we compare (c) useful 
squeezing and (d) tunnel coupling Q. Additionally we compare to results from the 
generic two- mode model, where a square- pulse is used (red dashed-dotted lines), (f) 
The two-body matrix elements of the OCT-solution are shown, similar as in figure [TT] 



also with figure [8|. In figure 12, panels (e) and (f), we report the tunnel coupling and 
two-body matrix elements. For the initial guess, which has a wildly oscillating density, 
also the tunnel coupling oscillates strongly and takes on a finite value after the control 
sequence. We interpret this as a signature of condensate excitations, which go beyond 
the two-mode MCTDHB model. In contrast, a smooth tunnel pulse and stationary final 
condensates are achieved for the optimal control. MCTDHB simulations with a higher 
number of modes indicate that the two-mode approach provides a very accurate level 
of description [63j. The two-body matrix elements of the optimized solutions show a 
complex behaviour and deviate from each other quite appreciably, which demonstrates 
that the dynamics depends critically on the orbitals [see equation (23) . 

For stronger interactions, we find that it becomes more difficult to decouple the 
condensate at the end of the tilting pulse, as shown in figure [13] for UqN = 1 and 
tphase = 5. Although wc casily achieve trapping of the orbitals in a stationary state, 
also for shorter Tp^ise oi Tphase, we find that number fluctuations are not constant after 
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Figure 13. Same as figure 12, but for U^N — 1, Tp^ise = 4 and Tphase = 5. (a) The 



inset magnifies ttie controls, (b)-(f) describe ttie same quantities as in figure [T2| 



the control sequence (results are not shown) , which is a signature of additional unwanted 
condensate excitations [63]. These excitations are most pronounced when the coherence 
factor tends to one, and the ungerade orbital becomes depopulated. In particular, we 
could not find pulses shorter than Tp^ise = 4 which lead to a final decoupling of the 



condensates. The same holds for the initial guess (blue dashed lines in figure 13). 
Although the condensate oscillations after the pulse sequence are moderate, they do not 
uncouple, as indicated by the growing tunnel coupling at later times. In contrast, the 
OCT control [green solid line in panel (a), magnified in the inset] achieves decoupling 
and a complete suppression of any tunnelling. Also the two-body matrix elements are 
stationary to a good degree. 

In conclusion, we find that optimized tilting pulses are crucial for avoiding 
condensate oscillations in the phase accumulation stage, and for reducing unwanted 
condensate oscillations. Adiabatic pulses might be orders of magnitude longer, and 
thus appreciably reduce the interferometer performance. 



5. Influence of temperature on the coherence of interferometer 
measurements 



Up to now we considered the atom cloud (BEC) at zero temperature. In realistic 
experiments the quantum system will be at some finite temperature, and especially for 
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the favorable low dimensional geometries temperature effects and decoherence due to 
thermal or even quantum fluctuations might become important [69l [TOl [TH [72]. We 
now turn to look at the effects of temperature on the decoherence of a split BEC 
interferometer. 



5.1. One- dimensional systems 

In one dimensional quantum systems (fc^T and ji < hujj_) fundamental quantum 
fluctuations prevent the establishment of phase coherence in an inflnitely long system 
even at zero temperature [70]. The coherence between two points in the longitudinal 
direction Zi and Z2 decays as \zi — ^2!""^^^^, where /C = nhy^ tlid / {qid'^) = ^"^idCh is 
the Luttinger parameter. Thereby C^^ = is the healing length, Uid is the 

ID density, and Qid = 2huj±as is the ID coupling strength in a system with transversal 
conflnement (jJ± and scattering length a^. For a weakly interacting ID system /C 1, 
and the length scale where quantum fluctuations start to destroy phase coherence is 

^,uant ^ _ (26) 

For weakly interacting one dimensional Bose gases (h is in the order of 0.1 to 1 fim. 
With /C > 10 one can safely neglect quantum fluctuations (see also [72j). 

These thermal phase fluctuations are also present in a very elongated 3D Bose 
gas. In such a flnite ID system one can achieve phase coherence (i.e., becomes 
larger than the longitudinal extension of the atomic cloud) if the temperature is below 
T^ = TiD^ ^ ^^loCh^z' To achieve a practically homogeneous phase along the BEC 
T < T$/10 is desirable [69] . 

For interferometry only the relative phase between the two interfering systems and 
its evolution are important. After a coherent splitting process, even a phase fluctuating 
condensate is split into two copies with a uniform relative phase. For interferometer 
measurements decoherence of this deflnite relative phase is adverse. 

The loss of coherence in ID systems due to thermal excitations was considered 
theoretically by Burkov, Lukin and Demler [73j, and by Mazets and Schmiedmayer [74j, 
and probed in an experiment by Hofferberth et al. [19j and Jo et al. [20j. The coherence 
in the system left is characterized by the coherence factor 



1 



^0 



(27) 



L 

where 0{z^t) is the relative phase between the two condensates. The angular brackets 
denote an ensemble average. The key feature in both calculations is that for ID systems 
the coherence factor ^^(/;) decays non-exponentially : 

^{t)ocexp[-{t/tof'] . (28) 
Burkov, Lukin and Demler [73] give for the characteristic time scale 



2.6l7r^ )C , (29) 
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whereas Mazets and Schmiedmayer [71] find: 



3,|if^V.3.2?!!#. (30) 



Even though the two show a different scahng with the Luttinger parameter /C, both are 
consistent with Hofferberth et al [19j within the experimental error bars in the probed 
range of /C ^ 30. 

There are two strategies to get long coherence times: 

• The characteristic timescale for decoherence scales like indicating that very 
long coherence times can be reached for low temperatures. 

• The time scale of decoherence scales with the ID-density as to oc nl^ in [74], 
and to oc n^(^ in f73]. Increasing the ID density will enhance the coherence time 
available for the measurement. 

Putting all together coherence times to in the order of lO^h/ fi can be achieved for T ^ 20 
nK and riio ^ 60 atoms//im. 

The above estimates have also to be taken with care, especially since the ID 
calculations are only strictly valid for T < /i, but we believe that they are a reasonable 
estimate also for elongated 3-d systems. 

The question is if the required low temperatures (T ^ /i/5 in the above example) 
can be reached in ID systems where two-body collisions are frozen out? In the weakly 
interacting regime, that is for small scattering length and relatively high density, this 
should be possible, because thermalization can be facilitated by virtual three body 
collisions [751 HHJ [77] . In the experiments by Hofferberth et al T ^ 30 nK and Uid ^ 60 
atoms / fim were achieved [22j 

Another point are dynamic excitations in the weakly confined direction. The 
interaction parameter Uq in general varies when A is changed and this can lead to 
longitudinal excitations. We find numerically that this effect is very small for the 
trapping potentials considered. In addition one can in principle account for the effects 
of the change in interaction energy by controlling the longitudinal confinement. 

5.2. Two-dimensional systems 

Quantum fluctuations are not important in 2D systems since a repulsive bosonic gas 
exhibits true condensate at T = 0. 

In a two dimensional Bose gas at finite temperature the phase fluctuations scale 
logarithmically with distance (r ^ Xt) [69j : 

(b(r)-^(0)])^^^ln(^^) , (31) 

where T2D = 27r/i^n2D/^ is the degeneracy temperature for a 2D system, n2D being the 
peak 2-d density. For T 3> /i the length scale is Xt ^ Ch^ where (h = ^/^/^^^zd^zd is the 
healing length in the 2D system with the effective coupling g2D- For weak interactions 



it is given by g2D = ^^f_^^^ with l± = \ h^^j^ [69]. For T <^ fi one finds the length scale 



CONTENTS 



29 



which is equal to the wavelength of thermal phonons. The phase coherence 
length is given by the distance where the mean square phase fluctuations become of 
the order of unity. 

h^Xr^wi^^ ■ (32) 

The decoherence in 2D systems was considered by Burkov, Lukin and Demler [73j , 
where they flnd a power law decay of the coherence factor: 

^{t) oc r^l^^KT ^ (33) 
for times larger then 

to - ' (34) 

where ji = n2Dg2D is the chemical potential of the 2D system. The temperature for the 
Kosterlitz-Thouless transition is given by 

TT n^-^ 1 

Trt o — 7^2D , (35) 

where n^^-^ < 7120 is the super fluid density. Again at sufliciently low temperatures the 
decoherence due to thermal phase fluctuations is smaller than the phase diflFusion due 
to the non-linear interactions during the phase accumulation stage. 



6. Summary and conclusions 

From our analysis of a double well interferometer for trapped atoms, it becomes 
evident that the main limiting factor to measurements with atom interferometers is 
the phase diffusion caused by the non linearity created by the atom- atom interactions. 
Consequently many of the recent experiments used interferometry to study the intriguing 
quantum many-body effects caused by interactions [181 UHl EHl 1221 [23] . 

Optimal control techniques can help improving interferometer performance 
signiflcantly by designing optimized splitting ramps and rephasing pulses, but the 
overall performance of the interferometer is still limited by atom interactions, and not 
by the readout, except for experiments with very small atom numbers. In general 
low dimensional conflnement of the trapped atomic cloud is better for interferometry. 
Nevertheless we found it difficult to get a performance for the minimal detectable shift 
AE'min < 10~^/x even for an optimized setting with a ID elongated trap. In addition we 
would like to point out that even though we did our analysis for a generic double well 
the same will hold for trapped atomic clocks [78] , where the signal comes from Ramsey 
interference of internal states. For the internal state interferometers the difference of 
the interaction energies is the relevant quantity to compare. 

The most direct way to achieve a much improved performance is to decrease the 
atom atom interaction. The best is to cancel it completely by either putting the atoms 
in an optical lattice, where on each site the maximal occupancy is 1, or by tuning 
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the scattering length = 0, which can in principle be achieved by employing Feshbach 
resonances [51] . Drastic reduction of phase diffusion when bringing the scattering length 
close to zero was recently demonstrated in two experiments in Innsbruck [52j and Firenze 
[53j. The big disadvantage thereby is that using Feshbach resonances requires specific 
atoms and specific atomic states. These states need to be tunable, and are therefore 
not the 'clock' states which are insensitive to external fields and disturbances. 

In an ideal interferometer one would like to use clock states, create strong squeezing 
during the splitting process by exploiting the non-linearity in the time evolution due to 
atom-atom interactions, and then, after the splitting turn off the interactions (by setting 
the scattering length to = 0). All together might be difficult or even impossible to 
achieve. 

For the interferometers considered here one can always reach low enough 
temperature to neglect decoherence due to thermal excitations even for ID and 2D 
systems. 

In addition to the interferometer scheme considered here, there exist other ideas 
of how atom interferometry could be improved. An interesting route will be to exploit 
in interferometry the correlations of the many boson states, and to establish a readout 
procedure which is immune to phase diffusion. One approach is to use Bayesian phase 
estimation schemes for the analysis of the phase sensitivity [791 ED] • Other proposed 
schemes include monitoring the coherence and revival dynamics of the condensates 
[81j . the measurement of a phase gradient along the double well potential containing 
tunnel-coupled condensates by means of a contrast resonance [82], or inhibition of 
phase diffusion by quantum noise [83j. For those more advanced ways to read out 
the interference patterns, however, the requirements for the temperature will become 
more stringent the better the new read-out schemes can compensate for the adverse 
effects of phase diffusion. 
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Appendix A. Optimality system for pulse optimization 

In this appendix, we give details about the optimality system for pulse optimization. For 
this purpose, we first discuss the pseudo-spin operators of equation ^ within MCTDHB. 
Since the many-boson state depends on both the number distribution and the orbitals, 
the operators now depend on position and time. Jz measures the difference between 
atoms in the left and right well, and Jx the difference between atoms in gerade and 
ungerade orbitals. Jy is then determined such that the pseudospin operators fulfill the 
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spin algebra. Then, 

Jx = -{a^gCig — a^ttg) , Jy = i{da)gae — d*a\ag) , = {da'gde + d*dldg) , (A.l) 

where d = J^^(j)g{x)(j)e{x) dx is the half sided overlap integral of the orbitals [33j. 
It measures the degree of orbital localization on the left-hand side. For our OCT 
optimization we define a cost functional for pulse optimization 

J(0„0e,C,A) = |^(aJ^^- (C(T)|J2|C(T)))' 

+ |[2 - mT)\<t>t)\' - \{MTM)\'] + I Ijm'dt, (A.2) 

where C is the atom number part of the wavefunction, and 71, 72, and 7 are weighting 
parameters. The first term quantifies the deviation of the actual AJy from the desired 
one, the second term accounts for trapping the orbitals in the ground state, and the 
third term penalizes rapidly varying controls and renders the problem well posed. In our 
OCT approach, the cost function is minimized subject to the constraint that the many- 
boson wavefunction fulfills the proper MCTDHB equations. Details of this approach 
can be found in Ref. [ 33j . 
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